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^^ ' We present a method which allows to include narrow-band correlation effects into the description 

of both valence and core states and we apply it to the prototypical case of nickel. The results of an 
ab-initio band calculation are used as input mean-field eigenstates for the calculation of self-energy 
"TX ' corrections and spectral functions according to a three-body scattering solution of a multi-orbital 

rr, , Hubbard hamiltonian. The calculated quasi-particle spectra show a remarkable agreement with 

photoemission data in terms of band width, exchange splitting, satellite energy position of valence 
states, spin polarization of both the main line and the satellite of the 3p core level. 
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PACS numbers: 71.10.-w, 79.60.-i 



tJ^ ■ I. INTRODUCTION 

-)-* 

^1 It is well established that the description of electronic states in narrow band materials requires improvements over 

G : the single particle approximation with a proper inclusion of on-site Coulomb interaction between localized electrons 
i-j^ lyi . In these systems the itinerant character of valence electrons which is clearly shown by the energy k-dispersion 

(— I I observed in photoemission spectroscopy coexists with strong local electronic correlation responsible of other observed 
Q ■ features such as satellite structures and band-narrowing effects. The interplay of localization and itineracy has also 
O ] been indicated as a possible explanation of the observed spin polarization of core level spectra through exchange 

coupling betweeen localized (core) and itinerant (valence) states [|[^ . 
^—{ [ In this paper we present a theoretical description of the valence and core electron states of nickel according to 

^ . a method recently developed which has been designed to treat highly correlated and highly hybridized systems |j] 
10 ' including both the itinerant character of band electrons and the strong localized electron-electron repulsion. This 
\l method allows to include narrow-band correlation effects into a first-principle band calculation; the single-particle 

^^ ' band states are determined according to the density functional theory in the local density approximation (LDA) and 
^^ \ the correlation effects are described as a 3-body scattering (3BS) solution of a multi-orbital Hubbard hamiltonian. 
r — This approach has been previously applied to the description of valence and conduction states of both model systems 
^>^ ' MM and realistic materials HRpI . We present here an extension of the method in order to treat both valence and 
core state on the same footing; this extension is possible in this scheme since it relies on a multi-orbital Hubbard 
C^ ' hamiltonian where core and valence states can coexist and on the 3BS method which can be applied for any value of 
effective on-site electron-electron repulsion. 

As far as the valence states are concerned various methods have been proposed to augment conventional band theory 
for the description of the electronic states of nickel; some of them are based on perturbative expansion either in the 



e-e interaction (GW approach [gj, second- order solution of the Hubbard hamiltonian [g^ 11 ij] ) or in the fluctuations 

of the electron occupation around LDA mean- field solution |1^ ; others apply a t-matrix scheme where the effect of 

electron correlation on one electron removal energies from a partially filled band is described as a hole-hole interaction 

^^ 0'Ell' This method is strictly valid only in the limit of an almost filled band (dilute limit) and its application to the 

S^ 1 case of nickel has been questioned |lO[ . In order to implement this approach it is necessary to include also electron- 

Vh " electron scattering channels and to solve a 3-body scattering problem involving two holes and one electron. This is 

5t 1 the spirit of the 3BS theory we apply here and which originally has been formulated by Igarashi [H|. This approach 

has recently been applied to the description of valence states of nickel 111] choosing, however, an approximate form 

of self-energy in terms of the antisymmetrized vertex function; in this way the self-energy turns out to be real, giving 

rise to peaks in the spectral density of unphysical zero width. Here we will instead adopt a version of the theory 

which avoids this shortcoming and which is based on the explicit solution of the 3-body scattering equations [|||3,0l . 

The interpretation of core level spectroscopies has up to now largely been based on atomic models which interpret 

the structures observed in the one electron removal spectra in terms of multiplet states formed by coupling the core 

hole to the unfilled valence shell [[l7|. T his scheme has in particular been applied to the photoemission spectra from 



the core states of transition metals [ 18[ attributing the observed characteristic line splitting to intra-atomic exchange 



interaction between core and valence electrons of an isolated atom. Such an approach has been seriously questioned 
since the observed splittings and the energy scale of the interaction are of the same order of magnitude as the valence 
band width ; a picture which takes into account the itinerant character of valence electrons seems therefore necessary. 
With the advent of spin polarized spectroscopies such as spin resolved x-ray photoemission y,|l9[| and magnetic circular 
dichroism EGlthe spin dependence of both the main line and the satellites of Ni core level spectra has been widely 
investigated [|19| , pi 22 2^ . We will show that the 3BS solution of a multi-band Hubbard hamiltonian, where full 



details of the valence band structure are included, can account for these spectroscopical features and interpret them 
in terms of on-site interaction between localized (core) and itinerant (valence) states. 

The paper is organized as follows: we present in sectio n pi| the multi-orbital Hubbard hamiltonian and its relationship 
to the band hamiltonian we want to implement; section |[lj describes the main characteristics of the method which we 
use to get an approximate solution of the Hubbard hamiltonian in terms of self-energy corrections to band eigenstates 
and of spectral densities; section |^ and ^ specialize to the case of valence and core states respectively; the results 



and the comparison with experiments are presented in section VI . 



II. MULTI-BAND HUBBARD HAMILTONIAN 

Band structure eigenvalues are in many cases good zero-order approximations to the excitation spectrum of a 
solid and it seems reasonable to use them as a starting point for the inclusion of correlation effects according to the 
Hubbard model; the implicit assumption is that among all the many body terms responsible of electron correlation, 
the Coulomb repulsion between electrons on the same site is the one which needs to be treated explicitly. To do this 
it is necessary to define precisely the relationship between band and Hubbard hamiltonians. Let us consider first a 
localized basis set (^iacr(r, s) with i labelling the localization site, a the orbital character, s and a the spin coordinate 
and eigenvalue respectively. The full many-body hamiltonian in second quantization is 

iatr a/3cr ij iajf3ljmS aa' 

with fiiacr = Ci^^Ciaa and Ciaa,Ci^„ dcstruction and creation operators. 

Here e^^^ and Uajp are the intra- and inter-atomic matrix elements of the one-particle hamiltonian ( kinetic energy 
+ ionic potential), while Viaa;jf3(T',i-ya-'.mSi7 are multi-center integrals involving the electron-electron interaction 

/2 
4'taai'^^ S)4'*jf3cr' (r', s')-r-^—yr<j>lja'{lc', s')(j)„rSa{r, s)drdr' 
ss- 

In this last expression the dominant contribution comes from the one-center integrals with i ^ j = I ~ m which are 
the usual on-site Coulomb term 

and exchange term 

The full many body hamiltonian can then be written as 

iaa a0a ij 



7 ,i^a0 ~ Jap) 2^ ntaaflifSa + 2_^ U^^p 2_^ ntaaflip-a 



2^ 

...(multi-center terms) (2-1) 



Different approximations of the exact hamiltonian (2.1) can be obtained using a mean field approach which amounts 
to neglect fluctuations in the electron occupation 

^iacT^i(3a' — ^iaa *^ ^ij3a' -^ '^ifJa' ^ ^iaa -^ "^ ^iacr -^ ^ ^ipa' -^ 
+ (n-iacr- < hiaa >) ' ("-i/3<T' ^ < "i/Scr' >) 
— fliaa < fliPcr' > +nipa' < fliacr > ~ < Uiaa >< fliPa' > 



where < .. > means a ground state average. The mean- field approximation can be applied to all the many body terms 
of {2A) transforming it into a single-particle hamiltonian 



H 



MF 



Ve^^-^n -I- 

/ , Haa "-laa- ^ 



apa ij 



(2.2) 



Any band structure calculation, where the interacting system is described as an effective single-particle problem, 
corresponds to the self-consistent solution of H^^^ . Another possibility is to apply the mean field approximation 
selectively to the multi-center integrals, keeping the full many body character in the one-center terms; in this way one 
gets a generalized Hubbard model 



a/3(T ij 



\Y. 



(2.3) 



ap Li <T i a 

Since H^^^ and H^ differ for the treatment of the on-site correlation - included in i/*^^ as a mean-field and treated 
as a many body term in H^ - it is easy to show that 

Due to the translational periodicity we can introduce an extended Bloch basis set 



(2.4) 



V', 



•l,{v,s) = -^^CV(k)e'^''-^'0.„.(r,. 



(2.5) 



and the corresponding relations for creation/destruction operators of electrons with wave vector k , spin a and band 
index n 



*k<T 



f-rj /_^^aa\'^)^ '' ^iaa 



at = ^Y.C:..{We-^^-^^cl„ 



Here C^(k(7) are the expansion coefficients of Bloch states in terms of localized orbitals, N the number of unit cells. 
H^ becomes 



kncT q/3 kk'p nn' mm' cr 

■ [u^f^C-A^yC^U^ + p)C^L.(k')*C^U(k' - p)aClaC+p.a™L,aI?lp_, 

+ (c/o/3 - Ja0)c:^{krc:'Ak + p)cj,(k')*<(k' - PKi^'^+^y^iK'-p. 



(2.6) 



Now e^no- includes also the kinetic part of the single-particle hamiltonian and the Coulomb and exchange integrals 
are assumed to be site-independent. In the same way H^^^ becomes 



H 



MF 






MF ;:,n^ -n 



ko-^ko- 



with 



MF _ H 



a/3 



_. occ _. occ 



k'n' 



k'n' 



(2.7) 

(2.8) 
(2.9) 



which is the analogue of eq. (2.4) for Bloch states. Notice that the sums over n' are over occupied states. Equations 
(2.8 2j| ) contain the correct recipe to include Hubbard correlation starting from band structure eigenvalues e^„^ and 
are essential in order to avoid double counting of e-e interaction. 



III. HOLE SPECTRAL FUNCTION, SELF-ENERGY AND THE FADDEEV METHOD 

We are interested in the hole spectral function 

^k.(^) = -^II^"i^"(k«^'^) (3-1) 

n 

which is the quantity directly related to the photoemission results we want to compare with. It describes the removal 
of one electron of wave- vector k, band index n and spin a and is related to the hole-propagator 



5-(kna,w) =-(*o|aklG(^)aCJ*o) ; z = -to + Eo{N,) +i5 



Eq (Ne) and |vl/o) define the ground state of the Ne particle system and 



(3.2) 



is the resolvent operator. The hole propagator can also be written in terms of the hole self-energy as 



g (kncr, w) - ^p 



^-ekna-^kn<T('^) 



(3.4) 



where Sr^^(a;) is the self-energy correction to band eigenvalues e^^^ . In order to calculate Sj^j^^(w) we proceed as 
in ref. [pj adopting a configuration-interaction scheme which consists in projecting the Hubbard hamiltonian on a set 
of states obtained by adding a finite number of e-h pairs to the Fermi sea, i.e. to the ground state j'l'o > of the 
single-particle hamiltonian. We will adopt the 3-body scattering (3BS) approach where this expansion is truncated 
to include just one e-h pair: the state whith one removed electron of momentum k and spin a is expanded in terms 
of the basis set including 1- hole configurations and 3-particle configurations (1 hole + 1 e-h pair) we will denote by 
\s) and |f) respectively 

\s) = akn<T |*o) 1^) = a]j3„3<^3aq2„20.2aqi„i^i |$o) (3-5) 

with 

qi + q2 - qa = k ; (Ti + 0-2 - era = cr 

To be consistent the basis set for the N-particle interacting system will include zero- and two-particle configurations. 
The ground state of the interacting iVg-particle system coincides then with the non-interacting one; this is obvious 
for the single-band hamiltonian discussed in ref. H (the state with zero and one e-h pair added coincide due to 
k-vector conservation) and still holds in the present case of multi-band hamiltonian since H has no off-diagonal 
matrix elements among 2-particle configurations. As a consequence in the 3BS approximation the hole propagator is 
the average of the resolvent over s tate s \s). 

By projecting the hamiltonian ( |2.6D over the complete set appropriate for the N^ — 1 particle system we get an 
approximate expression for H^ appropriate to describe one-electron removal 

HS^_^^Hi+H3 + V (3.6) 

Here Hi is associated to one-hole configurations 

Hi = {s\H"\s)\s){s\ 
H3 describes the contribution of 3-particle configurations 

w 
and V is the coupling between 1- and 3-particle states 



V = Y,{s\H"\t)\s}{t\ + h.c. 
t 

We leave the detailed expression of the matrix ele ments of the multi-band H^ to appendix A and proceed to sketch 



the method for the calculation of the resolvent (3^). We define the 3-particle resolvent, that is the resolvent associated 
to the 3-particle interaction 



and the Dyson equation which relates G{z) to it 

G{z)^F^{z) + F^{z)[Hi+V]G{z) (3.7) 

It is a matter of simple algebra to show that the hole propagator can then be expressed in terms of F3 as 

g~ (kna, io) = ~Gss{z) ^^^^ (3.8) 

w 

with the notation Gss = {s\ G \s), F^tt' = (^| -F3 1^') etc. Since the difference between the ground state energy of the 
A^e-particle system EQ{Ne) and the average of H^ over \s > states turns out to be 



Eo{Ne) - H,s - e^ncr + Qua - ^kri 



the band eigenvalues appear naturally in the denominator of the hole propagator and, comparing eq. (3.8) with ( p.4| ), 
we can identify the self-energy correction to band eigenvalues we are interested in as 

E- (kna, uj) = -J2 ^3tt' Vt>sV,t (3.9) 

w 

The determination of the hole propagator is then reduced to the calculation of F^tt' ■ This is done according to the 
Faddeev scattering theory |2J] as described for the single band case in reference [^ . The method consists in separating 
the 3-body hamiltonian in diagonal and non-diagonal parts 

iJ3^ = ^(t|iJ«|i)|i)(t| 
t 

H3^^ = 5](t|iJ«|0|i)(t'| 
W 

defining the diagonal 3-body resolvent 

pD (^) ^ — L^ 



and the scattering operator 



S = 4^^ + H^^'Fi'S 



The 3-body resolvent can then be written as 

F.^Ff+FfSFf (3.10) 

As shown in Appendix A and reference |6|] the non-diagonal 3-body interaction is the sum of two potentials 

^3 — ^h-h + Vh-e 

which describe h-h and h-e multiple scattering. 
We define partial scattering operators 



which are related to the scattering T-matrices 



through 



Sh-h - 


- Vh-h + Hfi^iiF^ S 


Sh-e = 


- Vh-e + Hh-eF.^ S 


itrices 

fh-h = 14- 


-h + Vh-hF^ Tfi-h 


fh-e = Vh- 


-e + Vh-eF^fh-e 


Sh-h = Th- 


-h + Th-hF^ Sh-e 



(3.11a) 



(3.11b) 



(3.12a) 



Sh-e — Th-e + Th-eFj, Sh-h 



(3.12b) 



These are the Faddeev equations whic h mus t be solved in order to get S — Sh-h + Sh-e, F^ from ( |3.1C| ) and 
Q~ (kna, w) from (3.8). Inserting (3.12) into ( 3.10| ) one gets the expression for the 3- particle resolvent in terms of 
scattering operators Sh-e and Th-h 



F.i 



Ff 



Th-h + Th-hF^ Sh-e + Sh- 



fF 



(3.13) 



Further steps are necessary in order to make this general method practical for real calculations - and further 
approximations as well. In the following we will specialize to valence and core states. 

IV. VALENCE STATES 

We adopt some semplifying assumptions for the description of photoemission from valence states of nickel. We will 
consider e-e correlation among d electrons only and neglect the orbital dependence of one-center integrals involving 
valence d electrons; we put 

Uap ~ Uvv for a, P — dorbitals 
— elsewhere 



and similarly for the exchange parameter JafS] moreover we will neglect Uyy 
the hamiltonian 2.t reduces to 



Jvv with respect to Uw In this way 



H 



H 



af3 kk'p nn' rnin' <y 

C^,(k)*C^:(k + p)C^_,(k')*C7r,(k' - p)a;^t 



EH -nt-n 
Ckn(T"ko-"kcr 

knCT 



k(T"k+pcr"k 



ra] 



-'k' — p— CT 



(4.1) 



Finally we will exclude configurations with e-h pair added to the majority-spin band as it would be strictly correct in 
the strong ferromagnetic limit where no empty states are available in the majority-spin band. 

The states which define the basis set for the A^e — 1 system are schematically illustrated in fig. 1; the non-zero 
interaction potentials responsible oi h — h and e — h scattering are also indicated. Notice that holes and electrons of 
parallel spin do no interact due to the assumption Uyy — Jyy ~ 0. 

We have already stressed that the Hubbard correlation enters the definition of quasiparticle energies twice: first as a 
mean- field correction to bare eigenvalues, transforming them into band ones according to eq. ( |2.8|) ; second through the 
addition of self-energy ( |3.9D . The calculation of this last quantity requires a generalization of the method illustrated 
in detail in ref. Q for the much simpler case of a single-band hamiltonian. The situation here is complicated by the 
sums over orbital indices appearing in the effective hamiltonian. This requires the definition of orbital-dependent 
diagonal Green functions and T-matrices as described in appendix B. As a result the self-energy correction to a band 
eigenvalue of wavevector k spin a and band index n turns out to depend on the quantum numbers n and k as 



Sk„.M = El^^V(k)l' 



T.^'^o^ 



empty 



^ |a^,(k')p-S0.(a.) 



(4.2) 



The k- vector and band-index dependence of self-energy is associated to the local orbital coefficients which modulate 
an orbital self-energy 

/"Ctt 

^pM -T. d' '^"-(^) ^^^hi^ - ^)[1 + UcpA'^^iu^ - e)] (4.3) 

where Haai^) is the spin-dependent orbital density of d valence states and T",".^ is the orbital dependent t- matrix 
describing the hole-hole multiple scattering 






Tf-ui^) - . , .;"".,, . (4.4) 



with 

gfioj) = r de' r & ""-^^)"^(^^l (4.5) 

A"^ includes the hole-electron scattering; it is determined by solving the integral equation 

/•oo 

A"P{uj, e) = B"^(w, e) + / de' n„-,(e') K''^{uj, e, e')A"^{u;, e') (4.6) 

JEf 



where 



K'^Piu;, e, e') = / de" n„_.(e") g',{uj + e" e)T^^^{uj + e")gUu; + e" e')T^^^{u: ~ e") (4.7) 



/ 



B"^(c., e)= / de' n„_,(e') g^{u; + e' - e)T,^^,(c. + e') 



X 



[gfiu; -e')+ de" n„_<,(e") g^,{^ + e' - e")gf{o. - e")r,"fjc. - e")] (4.8) 

JEf 

T^_g is the orbital dependent t-matrix describing the hole-electron scattering 

1 - Uapgi{uj) 
with 

na-o-(e)n/3o-(e') 



7-00 Jbj. w - e' ■ 



e — iJ 

- "J 

and finally 



(4.10) 



./_^ uj- e' -Id 



Equations (4.2 - 4.11) describe the procedure we have followed to calculate in practice self-energy corrections for the 



valence states of Ni reported in sec. VI, 



V. CORE STATES 

The localized character of core states is responsible for much stronger correlation effects associated with larger 
on-site Coulomb and exchange integrals. In order to adapt the hamiltonian ( |2.6D to the case of core states we make 
the following assumptions: we neglect the correlation among valence electrons and assume them to be described by a 



single band. In (2.6) the band index n can have just two values and as does the orbital index and we are left with a 



two-band haniiltonian 

H — Hyy + Hcc + Hcv (5-1) 

where H^^ is the band hamiltonian for valence electrons with the e-e interaction treated in the mean-field approxi- 
mation 



tlvv — / ^ ^kucr^kuQ-^k-uer 



kCT 



and 



ir _ Y^ H - 1 - I Y^ Y^ ^cc ., t ~ - 1 

rice - / ,ecCTQkca°kcCT ^ 2^ 2^ -T^ai^^^ak+pco-Ok/c-o-ak'-pc-o- 



N 

kcr kk'p (7 

Hcv = X! X! "/V \UcvCv-a(y)*Cv-a{\i' - p)aj^.^^ak+pcaaj^/^,^o-^k'-pD 
kk'p a 

+ yUcv ~ Jcv)^va{'^ ) ^va{'^ ~ PjOkco-'^k+pccrOk'iicr'^k' — pt)o 

Here Ucc = Jcc, Ucv, and Jcv describe the interactions between core - core and core - valence d electrons. 

Figure 2 shows schematically the configurations to be included for the 3BS description. Notice that configurations 
with one e-h pair added to the majority-spin band are now considered: even if the number of empty states available 
is small, the stronger value of the interactions makes these scattering channels no more negligible. As shown in fig. 2 
we have then to take into account scattering between particles of parallel spin proportional to Ucv — Jcv and between 
particles of opposite spin proportional to Ucv The extra configurations where all the particles (holes and electrons) 
have the same spin are called \z >. 

The procedure to calculate self-energy for core states is a straightforward extension of the one outlined for valence 
states. The core hole propagator turns out to be given by 

g- (kca, u;) = i — (5.2) 

uj - ef/ - Y, Fstt'Vt'sVst - J2 F3zz'V,',Vs, 

tt' zz' 

The presence here of extra configurations and extra interactions does not imply any major difference with respect to 



the case of valence states - just the addition of an extra term in the denominator of (5.2) and the necessity to solve 
separately two Faddeev problems to calculate F'nt' and F^zz' for opposite and parallel spin interactions respectively. 
This is a consequence of the fact that the hamiltonian ( |5.1[ ) does not mix \z > and \t > configurations. The Faddeev 
problem for the determination of F^zz' is solved in the same way as described in appendix B for F^tt' , substituting 
Ucv with Ucv — Jcv and parallel spin instead of opposite ones. Moreover the description of the core state in terms 
of a zero- width band with a delta function as orbital de nsity of s t ates drastically reduces the computational effort 



required to evaluate the non- interacting Green functions (4^, 4.10 , 4.11 ) 

The relationship between the energy e*^^ and the bare core eigenvalue involves as usual Q^^ (see eq. 2.8 and 2.9 
); explicitly one has 

ef5^ =eg+ [{Ucv - Jcv) < flva > +Ucv < nv-a >] + Ucc < flc-a > (5.3) 

It is interesting to make a comment concerning the role of the Coulomb repulsion Ucc- Since any 3-particle (1 hole 
plus one e-h pair) configuration must involve empty states, no multiple scattering is associated with Ucc and this 
quantity gives rise just to a mean field contribution. Core-valence interaction Ucv and Jcv on the contrary modify 
the bare core energies both through the mean field contribution and by self-energy corrections which originate from 
hole-hole and electron-hole scattering of opposite and parallel spin described by F^tt' and F^^z' respectively. 

VI. RESULTS 

The calculation of self-energy corrections requires as an input i) the mean-field eigenvalues and eigenvectors for 
valence electrons, ii)the energy of the core level and iii) the values of the Coulomb and exchange parameters Uw, 



Ucc^ Ucv, Jcv All these quantities can be deduced from an ab-initio LDA calculation; Coulomb integrals in particular 
can be obtained in the so-called constrained-density functional scheme ^^ - a procedure which is however not free 
of ambiguites (see for instance the discussion reported in ref. ||l^) and can lead to large variations in the estimated 
values. We have therefore adopted a mixed strategy, using results of ab-initio band calculations to get quantities i) 
only, treating all the others as free parameters. 

The band structure of ferromagnetic nickel has been calculated with the linear muffin-tin orbital (LMTO) method 
in the atomic-spheres approximation (ASA) including the combined correction term p6[. The tight-binding LMTO 



basis set |27 has been used, including 9 orbitals (s, p d) per atom. The resulting energy dispersion is shown in fig. 
3; the occupation numbers for valence d orbitals turns out to be < n^jj >= 4.68, < n„| >= 4.07. We have used the 
single-particle eigenvalues and the corresponding d contribution to eigenfunctions and orbital densities of states to 
calculate self-energy corrections and spectral functions according to the theory described in the previous sections. A 
value of Uyy — 2 eV has been chosen in order to reproduce the observed energy position of the valence band satellite; 
as we will show below we are able in this way to reproduce also other characteristics of the valence quasi-particle 
states such as band width, quasi-particle energy dispersion and exchange splitting; this is an important result and 
represents a success of the present approach: previous methods based on a semplified description of the scattering 



channels [14 15 1 in fact have not been able to reproduce at the same time the satellite energy position and the valence 
band width which turned out to be systematically overestimated for values of the Coulomb integral fixed to reproduce 
the satellite binding energy. 

Fig. 4 shows the comparison between our results and recent angle-integrated / spin-resolved photoemission data 
p|. We find that the calculated total spectral function D'^{uj) defined as 

kn 

closely reproduces the experimental energy distribution curves for each spin component; notice in particular that the 
6 eV satellite is observed clearly only in the photoemission from majority-spin states in agreement with our results. 
The two approximations we have adopted for the description of valence states, that is Uw — Jw — and the strong 
ferromagnetic limit make the self-energy corrections exactly zero for minority-spin bands and the comparison between 
our results and the spin-resolved experimental data confirms the validity of both these assumptions. As a further 
evidence of this we report in fig. 4 the minority-spin spectral function calculated after removing the assumption of 
strong ferromagnetism, i.e. considering also e-h pairs added to the majority-spin band; the small number of empty 
states and the relatively small value of the interaction Uyy make these scattering channels not efficient and the 
calculated spectrum is not significantly altered. 

It is possible to perform a more refined analysis by looking at the k- resolved spectral function; fig. 5 (a) shows 
the spectral function for the K point and e^^^ = -3.6 eV together with the corresponding self-energy E^^»(a;). The 
spectral function shows a quasi-particle peak plus a satellite: the first structure is associated to the pole of the 
hole propagator shown in fig. 5 (b) as the interception between ReCS^^^Juj)) and the straight line uj — e^^. The 

second structure is associated to the maximum of /?n(Sj^^t(w)) and as such it will occur at the same energy at any 
k-vector. This appears more clearly by extending the same analysis to the k-points along the high symmetry line of 
the Brillouin Zone and plotting the energy position of the maxima of D^^^{uj) = > ^knai^) ^'^ S*^^ ^^"^ quasi-particle 

n 

band structure of fig. 6. By comparing the quasi-particle band structure with the single-particle results it appears 
that the majority-spin eigenvalues are heavily affected by self-energy corrections, showing a strong reduction of the d 
- band width and the presence of the above mentioned 6 eV satellite. Since the majority-spin eigenvalues are shifted 
to lower binding energies while the minority-spin ones are unaffected, the splitting between majority and minority 
states becomes smaller than in the original single particle bands. This goes in the right direction since it is well 
known that single particle calculations overestimate this quantity; from fig. 6 it appears that the energy separation 
between majority and minority quasiparticle peacks around F point is 0.2 eV for the topmost band. This result is in 



remarkable agreement with a recent estimate reporting a split of 204±8 meV along the S direction |28|. The same 
is true for the quasi-particle energy dispersion as a whole, which is shown in fig. 7 compared with the results of 
angle-resolved spin-integrated photoemission spectroscopy of ref. ||29| . 

As far as core levels are concerned we consider here the 3p level of nickel as a test case. Since the core levels of an 
isolated atom are degenerate in spin their spin dependence is a solid state effect, associated to the interaction between 
the atom and the solid it is embedded in; in other words the spin dependence of core level energies is related to the 
spin polarization of the valence band on one side and to core-valence interactions on the other. According to our 
view the relevant quantities are the Coulomb and exchange integrals between core and valence orbitals which affect 
the bare atomic cor e en ergies both through the mean field term QJ^^ and self-energy corrections. The mean-field 



contribution of eq. (5_^) gives rise to a spin splitting proportional to Jcv and to the valence band spin polarization 



4]^ - ^c'l^ = Jcv (< nyi> - < fiy^- >) 
Notice that Ucv and Ucc do not enter this expression; they affect the bare core ievel e^ = e^, according to 

e^ + [Ucv (< fivT > + < fi^i >) + Ucc < fici >] = Cc (6.1) 

giving rise to a modified core energy level tc which is still spin-independent. The mean- field eigenvalues are related 
to this quantity as 

e^^^ - e - 7 <n > 

3/2 1 /2 

We have used spin-orbit split values of ec, that is Ec and tc , to fix the absolute value of the core level binding 
energy, choosing then Jcv to reproduce the spin-splitting of the main peak observed in the core-level photoemssion 
[g[. The last parameter, that is Ucv has been fixed in order to reproduce the satellites energy position. The values 
which optimize the agreement between our calculation and the experimental results are listed in Table | 

Fig. 8 shows the calculated self-energy for the creation of both majority and minority-spin hole in the core level 
3p^/^. The same analysis obviously applies to the other spin-orbit split level 3p^/^. In the case of the minority-spin core 
hole the self-energy presents two well defined structures: the one at lower binding energy is related to the scattering 
between particles of parallel spin with strength proportional to Ucv — Jcv (see fig. 2(f)) while the structure at higher 
binding energy is related to Ucv and to the scattering between opposite spin particles (see fig. 2 (e)). Notice that two 
independent factors determine the efficiency of the scattering process, the strength of the interaction on one side and 
the number of available states in the valence band on the other. In the case of a minority-spin core hole the weaker 
parallel spin interaction is compensated by the larger number of empty valence states and both the scattering channels 
involving parallel and opposite spin particles play a role. The same argument applies to the case of majority-spin core 
hole but now the weaker parallel spin interaction is associated to a small number of available empty valence states 
(see fig. 2 c); for this reason the interaction between opposite spin particles (see fig. 2 (b)) remains the only efficient 
scattering channel. As a consequence the self-energy for majority-spin core hole presents a single structure associated 

to Ucv 

As discussed in the previous section satellite structures are expected to occur at energies where the imaginary part 
of self-energy has a maximum giving rise to a complex pole of the hole propagator and therefore to a short lived 
excitation. The structures in the calculated self-energy we have just described and their origin are therefore essential 
in order to get a physical interpretation of the observed photoemission spectrum. We report in fig. 9 the calculated 
spin-integrated spectral density for the creation of a 3p core hole compared with the photoemission results of ref. [g2| . 
The spectrum shows a main peak at about 66 eV with two characteristic spin-orbit split structures (A, B) and two 
satellites (C, D). The decomposition of the spectral function into contributions from the two spin-orbit split levels 
3p^'^ and 3p^'^ and from different spins is also shown. The spectral function for the creation of a majority-spin 
core hole shows a main peack and a satellite for each spin-orbit split level; as discussed above in terms of the self- 
energy this satellite is associated to the only efficient scattering channel which comes into play after the removal of 
one majority- spin electron, that is to configurations of fig. 2. (b) where the majority-spin core hole interacts with 
opposite spin particles in the valence band, with a strength proportional to Ucv The spectral function for the creation 
of a minority-spin core hole presents instead two satellites for each spin-orbit split level: the one at the higher binding 
energy is associated to the configurations of fig. 2 (e) where the spin-down core hole interacts with opposite spin 
particles in the valence band with a strength proportianal to Ucv Due to the small number of empty states available 
in the spin-up band this satellites is less pronounced here than in the majority-spin spectrum. The satellite at lower 
binding energy is related to configurations of fig. 2 (f) and to scattering between parallel spin particles of strenght 

proportional to Ucv — Jcv 

The one-to-one correspondence between the configurations of fig. 2 and the satellite structures allows us to interpret 
them as as shake-up processes occurring after the removal of either a minority or a majority-spin core electron. The 
satellite at lower binding energy (C) is then associated to the creation of one minority-spin core-hole plus an e-h 
pair in the valence band of the same spin, giving rise to a bound state of three parallel spin particles which can be 
defined AS a triplet state; the satellite at higher binding energy (D) is related to the creation of either a majority or 
a minority-spin core hole plus a valence e-h pair of opposite spin giving rise to a singlet state. 

To analyze the spin polarization of the whole core hole spectrum in more details it is useful to consider the spin- 
resolved spectra and their difference D7 (uj) — D7 (w) reported in fig. 10. It appears that the spin polarity of structures 
A, B, C and D of fig. 9 is "down", "up", " down" , "up" respectively. This is again in agreement with what has been 



seen experimentally and reported in ref. [[19 22 12 
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VII. SUMMARY 

We have described a method for including short-range on-site interactions in the description of both valence and core 
states of a solid system. When applied to valence states of ferromagnetic nickel the method allows to get a quasiparticle 
band structure which compares much more favourably with the experimental observation than conventional mean- 
field LDA band structure eigenvalues , reproducing the observed band width, the energy dispersion, the satellite 
structure and the exchange splitting. Since the method does not rely on a perturbation expansion it has a wide 
range of application, including any correlation regime. The extension to core levels is quite natural and allows to 
take fully into account the itinerant character of valence electrons. We get then a physical picture of the 3p core 
level spectrum where the spin splitting of the main line is associated to the valence band spin polarization and to 
the core-valence exchange interaction; the two satellites are interpreted as arising from shake-up processes occurring 
after the removal of a core electron, involving the creation of a rf band e-h pair. Even though our present choice of 
empirically determining the parameters of the Hubbard hamiltonian ensures that we obtain an overall good agreement 
with experiments we believe that the possibility of reproducing both the satellite structures, the main line and their 
spin dependence with just four parameters can be seen as a non trivial result. The widely used atomic models - which 
use empirical parameters as well - require also an ad-hoc evaluation of the so called extra-atomic effects which allow 
to take into account the role of electrons on neighbouring atoms [£1[ ; such effects are on the contrary included here 
from the beginning since the full details of the valence band structure of the solid system are considered. In this sense 
the present approach can be considered an appropriate tool to describe the response of an itinerant electron system 
to the creation of a core hole. 



APPENDIX A: MATRIX ELEMENTS OF THE MULTI-ORBITAL HUBBARD HAMILTONIAN 



As discussed in section |^ the application of the 3BS me thod to the valence band states of nickel requires the 
explicit definition of matrix elements of the Hamiltonian 4.1 containing a non-interacting part that we call here iJ" 
and the on-site interaction among opposite spin electrons that we denote by H' . The matrix elements involve states 
\s > and \t > defined in 3.5 and pictorially depicted in fig. 1. We have 



k'n'a' 

{s\H°\t) = 

occ 
/i-\TTO\f\ _ \^ .H , H H H 

\i\ii \L/ — / ^ ^k'n'a' ^ ^qsns-a ^q2n2-(J qin^c 



k'n' a' 



\CUk)\ 



Q^ k"n" \k'n' 

{s\H'\t) = -E^C-,(g2)*Ci,(q3)C-(qi)*C^V(fc)4-,,.,.-,3 

a/3 
Q/3 l^ k'n' k" n" k'71' 

- KtAq2)\' E |<(fc')f - |c;;(9i)f E |^o-<.('^')f I 

k'n' k'n' ) 

(t'|Fli):=-E%^C'X(93)*C^::i.(93)C-(qi)*C';i(g;)5,,,,,<5„,,„-<5,,_,3,,;_,, 



+ 



a() 



aP 



E%^^a-(92)*C:L(92)C^;^('7i)*C;i(g;)<5,3,,,<5„3,„,<5,,+,„,- 



'l+92 



As discussed in section M the description of core states requires the inclusion of the interaction involving the 
exchange integral that we call here H" and the extension of the complete set to the three-particle configurations \z > 
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of fig. 2. It is easy to show that in this case the diagonal matrix elements are similar to the previous ones and that 
the non- zero off-diagonal ones are {s\H'\t) , {t\H'\t') , {z\H"\z') . 

APPENDIX B: FADDEEV EQUATIONS FOR THE MULTI-BAND HUBBARD HAMILTONIAN 



We illustrate the procedure which leads to eq. (p2) - ( 4.1C|) for the removal of one spin up electron. The relationship 



( p. 8]) between the diagonal hole-propagator and the 3-body resolvent is obtained by inserting the completness relation 

^|s><s|+^|i><t| 



into the identity 



z - H" 



that is 

[z - {H,)ss]G,s{z) - Y. ^stGtsiz) = 1 (Bl) 

t 

Since Gts{z) can be obtained by eq. ( ^ ) as 

Gtsiz)^J2^tt'Vr^,'^,s'iz) (B2) 



[z - {Hi),,]Gssiz) -J2y^tyt'.Fu'G,s +J2J2 V,tVt^,,Fu'G,s' = 1 (B3) 

tt' W s^s' 



one gets 



By substituting the explicit definition of the matrix elements of V it appears that the two summations in (|B3) involve 
terms of the kind 

with n = n' , n ^ n' in the first and second sum respectively. For this reason it seams reasonable to neglect in (|B3| ) 



the last summation with respect to the first one getting in this way eq. (3.8) as a result. 
According to eq. ( |3.9| ) the self-energy is given by the sum 



W 919293 nin2n3 Q/3 

■ E E E^"i(qi)^7i(q3)*crr(q2)c?T(k)* 

9i'J2?3 "i"2"3 f^ 

■ < qmi t (72^2 i qsns i WPi" + F^ {Su-e + Th-hF^ Sh-e)Fi' 

■ kWl t 'Z2'^2 i 93"3 i> 

where we have used the definition |3.5| for \t > and \t' >. Let us define now the orbital dependent free propagator 
describing h-h scatterings 



9ini92n2 ^ qiniT q2n2i q3n3i'' 



where the summation is over filled states of band indeces ni and n2. By using the definition ( 3.11 ) it is easy to show 
that 
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Yl C2]{mrC;i{q2r < qmi T q2n2 I 93^3 i \Gin-H\q[n[ t qW2 i qWs i> = 

1l92"2 



1 -f^Q/333 (qsw) 

It is also useful to define the orbital-dependent t-matrix for h-h scatterings 

Th-M-^n^^) = ^ 

Similar definitions and relations hold for e-h scatterings. It is then a matter of simple algebra to transform eq. (B5) 
into the form 



^i^3U'^^t'aKt=^|C^'(kT)|' 



a qana 



where 



gi92 nin2 

• E E <Hqi)<^(q3)*^rr(q2) 



■ < qin-i T <?2"2 i g3'^3 i \FzSh^eF^\q[n'-^ t '72'^2 i 'J's^-s ■!> 

From now on the procedure is the same as the one described in ref. H, adopting in particular the so-called local 
approximation |IQ] which allows to transform all the summation over k-vectors into integrals involving the density of 
states. 
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FIG. 1. Schematic representation of the basis set for the configuration expansion of the interacting state with one electron 
removed from the majority-spin band. V , Vh-h, Vh-e describe the coupling between 1- and 3-particle states, the hole-hole and 
the hole- electron interaction respectively. 

FIG. 2. Schematic representation of the basis set for the configuration expansion of core states. 

FIG. 3. Single particle band structure of nickel obtained using the LMTO method. Energies are referred to the Fermi energy. 
Open triangles, minority-spin states; filled triangles, majority-spin states. 

FIG. 4. Density of quasi-particle states of nickel for majority-spin (a) and minority-spin (b) spin compared with experimental 
results of angle-integrated spin-resolved photoemission results (filled triangles) of ref. S . The results for minority-spin bands 
obtained without the strong- ferromagnetic-limit approximation is shown as a broken line in panel (b). 

FIG. 5. Spectral function (a) and self-energy (b) for the K point and e^J}^ = -3.6 eV. The interception between 
Re (E~(kn Ti'^)) a-nd the straight line ui — e^^^ indicate the position of the quasi-particle pole. 

FIG. 6. Quasi-particle band structure of nickel along the high symmetry directions of the Brillouin Zone. Energies are 
referred to the Fermi energy. Open triangles, minority-spin states; filled triangles, majority-spin states. 

FIG. 7. Comparison between the calculated dispersion of hole quasi-particle states (circles) for majority-spin bands and 
angle- resolved spin-integrated photoemission results (diamonds) of reference [E9| . 

FIG. 8. Self-energy for the creation of majority-spin (a) and minority-spin hole (b) in the core level 3p^'^. 

FIG. 9. Spin-integrated spectral function for the creation of a 3p core hole compared with the photoemission results of ref. 
[h| (a) and its decomposition into contributions from the spin-orbit split level Zp^'^ (continuous line) and 3p''" (dashed line). 
Panel (b) refers to majority-spin state and panel (c) to the minority-spin one. 

FIG. 10. (a) Calculated spectrum for the creation of a majority-spin (continuous line) and minority-spin (dashed line) 3p 
core hole, (b) Difference spectrum between the two spin component. 

TABLE I. Parameters used in the calculations of Ni 3p core hole spectrum 

(--'Ctl 'J CV Cc Cc 

5.00 eV 2.50 eV -59.00 eV -57.50 eV 
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